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Selfsynchronization and dissipation-induced threshold in collective atomic recoil lasing 
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Networks of globally coupled oscillators exhibit phase transitions from incoherent to coherent 
states. Atoms interacting with the counterpropagating modes of a unidirectionally pumped high- 
finesse ring cavity form such a globally coupled network. The coupling mechanism is provided by 
collective atomic recoil lasing (CARL), i.e. cooperative Bragg scattering of laser light at an atomic 
density grating, which is self-induced by the laser light. Under the rule of an additional friction 
force, the atomic ensemble is expected to undergo a phase transition to a state of synchronized 
atomic motion. We present the experimental investigation of this phase transition by studying the 
threshold behavior of the CARL process. 
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Lasing mediated by collective atomic recoil (CARL) 
has been predicted as the analogue to free-electron las- 
ing Q and observed recently 0. In the experiment an 
ensemble of cold atoms couples to the two counterprop- 
agating modes of a unidirectionally pumped high-finesse 
ring-cavity and collectively scatters photons between the 
modes. A collective instability leads to self- amplification 
and to exponential gain for a light mode on one hand 
and atomic bunching on the other. In the absence of en- 
ergy dissipation for the kinetic degrees of freedom, the 
gain must remain transient. The introduction of a fric- 
tion force for the kinetic energy of the atoms permits 
however a steady-state operation of the CARL at a self- 
determined frequency (viscous CARL) 0. 

The CARL represents a system of coupled oscillators. 
Collective dynamics in networks of weakly coupled sys- 
tems is a very general phenomenon. The richness of this 
paradigm is illustrated by the large amount of examples, 
ranging from physical systems, like arrays of Josephson 
junctions or phase-locked lasers, to biological systems like 
cardiac pacemaker cells or chorusing crickets 0. Ku- 
ramoto, Strogatz and others 0, have shown that en- 
sembles of coupled oscillators operating at different fre- 
quencies or being subject to stochastic noise synchro- 
nize, if their number and their mutual coupling strength 
exceeds a critical value. I.e. cooperative action starts 
beyond a threshold value of the coupling strength after 
crossing a thermodynamical phase transition. 

In the present paper, we investigate the threshold be- 
havior of the CARL as an example of a phase transi- 
tion of the Kuramoto type 00: If, for viscous CARL, 
the dissipation (or cooling) mechanism is limited to finite 
temperatures by some diffusion process, a phase transi- 
tion occurs in the atomic density distribution when the 
system is pumped at threshold. Consequently, a min- 



imum pump power is necessary to start CARL lasing, 
just like in ordinary lasers. However, in extension to 
the original Kuramoto model, the collective frequency is 
self-determined. We report on the observation of such a 
threshold and characterize it in terms of the theoretical 
model presented in Ref. 0. 

The optical layout of our experiment has been out- 
lined in Ref. 0. A titanium-sapphire laser is tightly 
phase-locked to one of the two counterpropagating modes 
of a ring cavity by means of a Pound-Drever-Hall type 
servo control. The amplitude decay rate of the ring cav- 
ity is k — (27r)22 kHz. In the following, we label the 
modes by their complex field amplitudes scaled to the 
field per photon a±. The intracavity light power is then 
P± = hui5\a±\ 2 , where 6 — 3.5 GHz denotes the free 
spectral range of the cavity. The phase dynamics of the 
two counterpropagating cavity modes is monitored via 
the beat signal between the two outcoupled beams: Any 
displacement of the standing wave inside the ring cavity 
is translated into an amplitude variation of the observed 
interference signal, Pbeat = ftwS\a+ + a_| 2 . 

The contrast of the standing wave Pbeat is weak as 
compared to the pump power (a few %). Assuming a+ 
real and |cv_ | <C a+, the probe beam power can be related 
to the contrast of the beat signal AP cont = 4fiw<5a + |a_|, 



P_ = 



AP 2 



16P+ 



(1) 



In our experiment, we load 85 Rb atoms from a stan- 
dard magneto-optical trap (MOT) into the optical dipole 
potential, which is generated by a TEMoo mode of the 
unidirectionally pumped ring cavity and tuned to the red 
of the D\ line. Typically 2 x 10 6 atoms are trapped and 
form a basically homogeneous 4 mm long cloud along 
the cavity axis around the waist of the cavity field. The 
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cloud reaches a peak density of about 2 x 10 9 cm" 3 and 
a temperature of a few 100 /zK. Absorption images of the 
atomic density distribution taken after a time-of-flight 
and spectra of the velocity distribution obtained by excit- 
ing recoil-induced resonances (MR) |8( yield information 
on the kinetic degrees of freedom of our system, which is 
complementary to that on the field amplitudes a±. 

A genuine problem of the CARL is the following: 
In the absence of damping for the external degrees of 
freedom the CARL process continuously accelerates the 
atomic center-of-mass 0, ■ Even though the acceler- 
ation decreases because the Doppler-shifted CARL fre- 
quency eventually drops out of the cavity resonance, it 
never reaches a stationary value. In fact, being focusscd 
on studies of transient phenomena, the original CARL 
model Q does not consider relaxation of the translational 
degrees of freedom. On the other hand, standard meth- 
ods of optically cooling atoms are based on controlled 
dissipation, e.g. optical molasses. Close to resonance 
the motion of atoms in an optical molasses is well de- 
scribed by a friction force. In our experiment, we har- 
ness this dissipation mechanism and subject the dipole- 
trapped atomic cloud to an optical molasses. We use the 
laser beams of the MOT and tune them 50 MHz below 
the cooling transition (D-2, F = 3 — > F' = 4). In this 
situation, the beat frequency oscillations quickly reach a 
stable equilibrium frequency between Auj/2tt = 100 kHz 
and 170 kHz, which corresponds to an atomic velocity of 
7 to 13 cm/s. 

In order to observe a threshold behavior in experi- 
ment, we adiabatically ramp up and down the inten- 
sity of the pump laser. The CARL radiation is mon- 
itored by recording the time evolution of the beat fre- 
quency of the counterpropagating modes. The curves 
shown in Fig. ^a) represent frequency spectra obtained 
by Fourier transforming the beat signal restricted to suc- 
cessive time-intervals. The peaks' locations then denote 
the instantaneous frequency shift of the probe beam, and 
their heights reflect the standing wave's contrast. The 
pronounced dependence of the CARL frequency on the 
pump intensity revealed by the series of Fourier spectra 
is emphasized in Fig. QJb) ■ The intensity of the CARL 
radiation shown in Fig.^c) decreases with the pump in- 
tensity, but more important is the fact that the curve 
exhibits a minimum pump intensity required to initiate 
CARL lasing. Just like in common lasers, only if the en- 
ergy fed to the system exceeds the losses, the laser emits 
coherent radiation. 

The experimental observations may be discussed in 
various ways. The model of Ref. 2j to describe the 
impact of optical molasses on CARL consisted of sim- 
ply adding a friction force, proportional to a coefficient 
jf r , to the equations governing the atomic dynamics. 
This procedure certainly represents a coarse simplifica- 
tion. For example it predicts that the atoms quickly 
bunch under the influence of the molasses, and are cooled 
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FIG. 1: (a) Sectionwise Fourier-transform of the interfer- 
ence signal Pbeat with the pump laser power being ramped 
up and down. The dotted line is proportional to the pump 
laser power, (b) Dependence of the CARL frequency on the 
intracavity pump power, (c) Dependence of the probe field 
intensity on the pump power. The CARL laser threshold is 
around P+ — 4 W intracavity power. The fitted curves are 
based on a Fokker-Planck theory outlined in 0. (d) Cal- 
culated bunching parameter. The parameters are jf r = 4k, 
N = 10 6 , A a = -(2tt)1.7 THz, and T = 200 \iK. 



until the temperature of the cloud is T = 0, and it de- 
nies the presence of any threshold. In reality, the mo- 
lasses temperature is limited by diffusion in momentum 
space, i.e. heating. To account for this heating, one may 
supplement the dynamic equations for the trajectories of 
individual atoms (Ref. |2(, Eq. (1)) with a stochastically 
fluctuating Langevin force £ n (t) with (£ n (t)) = and 
(£n(i)£m( T )) — 2lf r D5 mn 5(t— t), where the diffusion co- 
efficient D — a 1 1 'jf r is proportional to the atoms' equi- 
librium temperature, which is related to the Maxwell- 
Gaussian velocity spread by a = 2kyJksT /m: 

n = 4eiU a+ (a_e~ ie " - a*_e l9n ) - 7/r#n + &, • (2) 

Here we defined 9 n = 2kx n as the position of the n th 
atom along the optical axis normalized to the optical 
wavelength and assumed the pump laser stabilized on 
resonance with the cavity, a + is set real and constant. 
N is the total atom number, and s = Hk 2 /m is twice 
the recoil frequency shift. The coupling strength U a (or 
single-photon light shift) is related to the one-photon 
Rabi-frequency g and to the laser detuning from reso- 
nance by Uo = g 2 /A a - The functional dependence of 
the quantity a_ on the order (or bunching) parameter 
b = \b\e 1 ^ = N^ 1 J2 m e%6m i s determined by the differen- 
tial equation 



-Ka_ — iNUoo:+ b . 



(3) 



An alternative to simulating trajectories of individual 
atoms is to calculate the dynamics of distribution func- 
tions. Particularly adequate to the problem of diffusion 
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in momentum space induced by optical molasses is a 
Fokker-Planck approach. Here the thermalization of the 
atomic density distribution P towards an equilibrium be- 
tween cooling and heating is described by the interplay 
of friction and diffusion pj. In the limit of strong vis- 
cous damping, where we may adiabatically eliminate the 
atomic momenta by setting 9 n = 0, the Fokker-Planck 
equation associated to the Langevin equation (J2J) reads: 



dP _ 4ieU a+ d (a*_e 2i0 - a_e 
~dt ~ 



P 



1th 



06 



d 2 P 
•^,(4) 



and the bunching parameter is given by b — Pe lS d9. 
The solid lines fitted to the data in Fig. ^ are calculated 
by numerical integration of this Fokker-Planck equation. 
Fig. njd) shows the calculated evolution of the bunching. 
Apparently, the bunching vanishes below the CARL las- 
ing threshold and tends towards 1 as the pump power is 
increased. The threshold behavior of the radiation mode 
is thus intrinsically connected to atomic self-organization. 
As has already been realized by Kuramoto, the Fokker- 
Planck equation predicts the occurrence of a thermody- 
namical phase transition. An alternative approach devel- 
oped by Bonifacio and Verkerk jlfj studies the evolution 
of the atomic phase-space distribution described by the 
Vlasov equation towards equilibrium with a single rate 
7/ r , and Javaloyes et al. |llj found that the Vlasov ap- 
proach leads to a phase transition. 

Ref. [7j pointed out that the threshold should depend 
on various parameters like the atom number, the cou- 
pling strength, the friction coefficient and also on the 
atomic temperature. However a proper scaling ,73 shows 
that the threshold is ruled by only two independent quan- 
tities. Analytic expressions for the threshold conditions 
are obtained from a linear stability analysis of the Fokker- 
Planck equations. As shown in Ref. Q, by introducing 
the CARL parameter p = (NU^a 2 h /2e 2 ) 1/3 , the condi- 
tion for lasing is given by 



KG 



7/r 



(<j 2 + K j fr ) 2 < (2e P f 



(5) 



In the good cavity limit K^f r < er 2 , we may neglect the 
second term in the bracket. At threshold, where the 
equality in the above equation holds, the CARL param- 
eter becomes p t hr = {i^/lfr) 1 ^ (cr/2s). Therefore the 
threshold pump power Pthr - 



l +,thr 



7/r 



FihjSa 2 , thr follows with 

V2 „3 



(6) 



The threshold power and the threshold CARL frequency 
are related by [7j, Aw^ r = {2ep t hrf' ft/7/ r - As a conse- 
quence, the CARL frequency is independent of the atom- 
field coupling, 



Auj 



thr 



7/r 



1/2 



(7) 
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FIG. 2: Pump power threshold measured for various laser 
detunings (a) and atom numbers (b) determined from curves 
like those shown in Fig. Q (c) CARL frequency shifts as a 
function of the laser detuning, (d) Pump power threshold 
plotted as a function of the temperature derived from the 
threshold CARL frequency according to Eq. 0. 



Figs.|2a) and (b) show measurements of the threshold 
pump power as a function of the laser detuning and the 
total atom number. Both experimental curves have been 
obtained from the same data set P + , t h r (N, A a ) scaled 
either to an arbitrarily chosen reference atom number 
N = 10 6 (curve (a)) or to a reference laser detuning 
A a = — (27r)1.5 THz (curve (b)) using the relationship 
©. The fitted curves are obtained from the same equa- 
tion @ by adjusting the temperature to T — 200 pK 
and the friction coefficient to 7/ r « 4k. Fig.^c) demon- 
strates that the CARL frequencies measured at thresh- 
old do apparently not depend on the coupling strength. 
The horizontal line indicates the frequency correspond- 
ing to the temperature T — 200 pK according to Eq. J7J . 
On the other hand, theory predicts a dependence of 
the threshold power on the atomic temperature: Hot- 
ter atomic clouds have a broader velocity distribution a 
and exhibit a higher threshold power. Even though the 
temperature data shown in Fig.EId) seem to confirm this 
trend, they are too uncertain to be used to improve the 
fits of Figs.EIa) and (b). 

The steady-state operation of our system is ensured 
by the optical molasses. As soon as the molasses is 
turned off, the equilibrium is lost and the atoms and the 
standing-wave accelerate each other, provided the atoms 
are bunched. Fig.|^a) demonstrates the acceleration pro- 
cess. The instantaneous CARL frequency and intensity 
are again obtained from a sequence of Fourier spectra 
taken over successive time periods. The evolution of the 
probe power calculated with Eq. Q is shown in Fig.[3Jb). 
At times t < 0.5 ms, when the molasses is present, the 
CARL frequency is fixed at 170 kHz. As soon as the 
molasses is switched off, the probe's amplitude dimin- 
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ishes while its frequency detuning from the cavity res- 
onance increases. The acceleration indicates the occur- 
rence of bunching induced by the optical molasses. The 
self-organization of the atomic density distribution out 
of a homogeneous cloud into a spatially ordered arrange- 
ment spontaneously breaks translational symmetry and 
represents the strongest signature of CARL action Q . 
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FIG. 3: (a) Evolution of the CARL frequency after molasses 
has been switched off at time t = 0.5 ms. (b) CARL power 
calculated via Eq. JTJ. The fits are based on the theoretical 
formulae Eqs. (2) and (3) of Ref. 0. 

The probe field a_ depends on the location of all 
atoms in a collective way. Consequently, Eq. J5J) de- 
scribes a mean-field type dependence of the location of 
every single atom on all the others. This fact gets par- 
ticularly transparent in the limit of strong viscous damp- 
ing, 9 n — 0, if we furthermore assume that at steady 
state the optical standing wave propagates with a con- 
stant amplitude at a constant velocity. This condition is 
formulated by d_ = iuj ca a^, where the probe beam fre- 
quency shift uj c(1 = io ca (b) depends on the atomic bunch- 
ing. With these approximations we substitute the solu- 
tion of Eq. J2J into Eq. © and obtain, defining the cou- 
pling constant K = %eNU^a\u!ca/ (lfr{^ca + k2 )) an d 
restricting ourselves to the good cavity limit k <C uj ca . 



7/r 



+ K\b\sm {ip-0 n ) 



(8) 



The equation © describes the dynamics of N coupled 
oscillators with fictitious frequencies uj n = Cn/jfr- This 
system, which has been investigated by Kuramoto Q, 
predicts the synchronization of those oscillators over 
time, whose frequencies satisfy uj n < K\b\. The anal- 
ogy between our viscous CARL and an ensemble of self- 
synchronizing harmonic oscillators resides in the follow- 
ing correspondences: The phases of the oscillators are 
represented by the positions of atoms. Synchronization 
of the oscillators corresponds to bunching of the atoms. 
The role of friction is to provide a steady atomic center- 
of-mass velocity to which the individual atomic velocities 
may lock. In the case of CARL and unlike for the Ku- 
ramoto model the collective oscillation frequency is self- 
determined. Diffusion is the source of disorder, which 
rules the phase transition by competing with the dynam- 
ical coupling, in contrast to the Kuramoto model, where 
disorder occurs via distributed natural frequencies. 

To conclude we point out, that the viscous CARL sys- 
tem is representative for a vast class of systems. Under 



the rule of optical molasses the coupled field-atom sys- 
tem constitutes an ideal model system for an ensemble 
of weakly coupled oscillators. Despite the fact that the 
details of the molasses dynamics are complicated, its im- 
pact on our CARL system is fairly well described by two 
constants, the friction and the diffusion coefficient. Al- 
though the system is purely classical, we deal with micro- 
scopic particles. This bears the possibility of transferring 
the system to the quantum regime, and thus to study 
the coupling of large ensembles of quantum oscillators. 
Furthermore, classical networks of dynamical systems in 
general depend much on the details of how the coupling is 
realized. In contrast, the coupling in our system is gener- 
ated by the fundamental interaction between atoms and 
light, which is very well understood and even controllable 
by experiment, e.g. via the tunable friction force or the 
laser detuning. Because it is mediated by a delocalized 
object, i.e. a standing light wave, the coupling is instan- 
taneous (no retardation effects) and truly uniform (every 
atom couples with the same strength to all its neighbors). 
The large number of oscillators ensures that the system 
is in the thermodynamic limit. 
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